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We study the transport of electrons through a long quantum wire connecting two bulk leads. As 
the electron density in the wire is lowered, the Coulomb interactions lead to short-range crystalline 
ordering of electrons. In this Wigner crystal state the spins of electrons form an antiferromagnetic 
Heisenberg spin chain with exponentially small exchange coupling J. Inhomogeneity of the electron 
density due to the coupling of the wire to the leads results in violation of spin-charge separation 
in the device. As a result the spins affect the conductance of the wire. At zero temperature the 
low-energy spin excitations propagate freely through the wire, and its conductance remains 2e'^/h. 
Since the energy of the elementary excitations in the spin chain (spinous) cannot exceed tvJ/2, the 
conductance of the wire acquires an exponentially small negative correction SG oc — exp(— 7rJ/2r) 
at low temperatures T <^ J. At higher temperatures, T ^ J, most of the spin excitations in the 
leads are reflected by the wire, and the conductance levels off at a new universal value e^ /h. 

PACS numbers: 73.63.Nm, 73.21. Hb, 75.10.Pq 



I. INTRODUCTION 

The quantization of conductance of one-dimensional 
(ID) electron systems in units of 2e^ /h was first observed 
in experiments with quantum point contactsi^ The lat- 
ter consist of a short (well under 1 fj,m) ID constric- 
tion connecting two bulk two-dimensional leads. Further 
progress in fabrication of low-disorder devices resulted 
in observation^ of similar conductance quantization in 
quantum wires of several microns in length. Experimen- 
tally the quantization is observed as very flat plateaus in 
the dependence of linear conductance on the voltage at 
the gate controlling the electron density in the wire. 

In a number of recent experiment s^'^i^i'^i^i^i-'^^'-'^-'^i"'^^'"'^? de- 
viations of conductance from perfect quantization have 
been observed at low electron density. These devi- 
ations manifest themselves as negative corrections to 
the conductance at the beginning of the flrst quantized 
plateau. The correction is usually small at the low- 
est temperatures available, but becomes significant at 
T ~ IK. In typical samples^*Sii the conductance levels 
off at high temperatures and forms a quasi-plateau at 
about 0.7 X (2e^//i). This phenomenon is often referred 
to as the 0.7 structure. Despite the numerous theoret- 
ical attemptai4.i-^ii6.i^ii«.i9,2o,2i,22^23,24,25,26 at the inter- 

pretation of the 0.7 structure, its origin remains unclear. 

The analysis of the existing experimental data shows 
that the 0.7 structure is sensitive to the length of the 
one-dimensional region connecting the two-dimensional 
leads. The structure tends to be relatively weak in the 
short contacts?^ where no quasi-plateau is observed even 
at high temperatures. On the other hand, in longer 
samplesifliiiiiSii^ the plateau is observed even at the low- 
est temperatures available. The effect is also somewhat 
stronger, with the quasi-plateau moving to a lower value 
of conductance G ~ 0.5 x {2e^/h). 



In this paper we consider conductance of a long quan- 
tum wire in the regime of low electron density n. Focus- 
ing on the case of a GaAs device, we assume quadratic 
energy spectrum for free electrons e(p) — p'^ /2m, where 
m is the effective mass of the electrons. The typical ki- 
netic energy of an electron at zero temperature is of the 
order of the Fermi energy, Ep = {■kTiti)^ j'&m. It is im- 
portant to note that at low density n <g; a^^ the kinetic 
energy is small compared to the typical energy e^n/e of 
Coulomb interaction between electrons, where e is the 
dielectric constant, and as — eh^ /me^ is the effective 
Bohr radius of the material. Thus in the limit of low 
density the electrons can be viewed as classical particles 
placed at equidistant positions to minimize the Coulomb 
repulsion. Such a picture was first proposed by Wigner— 
and will be referred to as the Wigner crystal. 

Although the quantum fluctuations of electrons near 
their equilibrium positions destroy the long-range order 
in the Wigner crystal, its short-range structure strongly 
affects the transport through the wire. In particular, the 
electrons occupying well-defined sites of a Wigner lattice 
can be viewed as an antiferromagnetic spin chain with 
exponentially small exchange constant J. The appear- 
ance of a new energy scale J <C Ep significantly affects 
the physics of the electronic transport through the wire. 
This effect is most important at intermediate tempera- 
tures, J <C T <C Ep, where it results in a considerable 
suppression of the conductance of the wire. 

The physics of this phenomenon is controlled by the 
effect of spin-charge separation in one-dimensional inter- 
acting electron systems. The latter refers to the fact that 
fermionic quasiparticles cannot be viewed as elementary 
excitations of the system, i.e., it no longer behaves as a 
Fermi liquid. Instead, the system displays Luttinger liq- 
uid behavior, and the elementary excitations are bosonic 
waves of charge and spin densities propagating at dif- 
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ferent velocities. As a result, when an electron enters 
the one-dimensional region from one of the leads, it is 
decomposed into charge and spin waves. At low temper- 
ature T <^ J both waves pass through the wire, and upon 
reaching the other lead they reassemble into a fermionic 
quasiparticle. This process can be interpreted as perfect 
transmission of electrons through the wire and gives the 
standard value of conductance G = 2e'^/h. On the other 
hand, at T 3> J, the bandwidth ~ J of the spin ex- 
citations in the wire is small compared to their typical 
energy ^ T. As a result, only the charge excitations pass 
through the wire, whereas the spin ones are reflected back 
to the lead. We show below that this additional scatter- 
ing of spin excitations by the wire reduces its conductance 
to e'^/h 

In Sec. ^ wc study the applicability of the Luttinger 
liquid description to the ID Wigner crystal and show that 
it is only valid at energy scales below J. On the other 
hand, the property of spin-charge separation is more gen- 
eral and persists at energy scales above J. We review the 
known results for the conductance of a quantum wire in 
Sec. mil The most important consequence of the spin- 
charge separation in quantum wires is the conclusion that 
spin degrees of freedom do not affect the conductance, 
which remains quantized at 2e'^/h. In Sec. II VI we show 
that the spin-charge separation is violated when the wire 
is connected to the leads. As a result the spin subsys- 
tem affects the propagation of electric charge through the 
wire and contributes to its resistance. This contribution 
is studied in Sec. El where we find that at J <C T the con- 
ductance of the device reduces from 2e^ /h to /h. The 
relation of our results to experimental measurements of 
conductance of quantum wires is discusses in Sec. IVII A 
brief summary of some of our results has been reported 
in Ref.|2g 



II. SPIN-CHARGE SEPARATION IN 
QUANTUM WIRES 

It has been known since the seventiesSS. that the low- 
energy excitations of a ID system of interacting electrons 
are the charge and spin waves propagating independently 
of each other at different velocities. This result is valid 
at the energy scales low compared to the bandwidths Dp 
and of the charge and spin excitations. In the case of 
not very strong interactions both bandwidths are of the 
order of the Fermi energy, and the picture of completely 
separated charge and spin excitations is appropriate at 

At low electron density n ^ the interactions be- 
tween electrons are strong. Wc show below that as a re- 
sult the velocity of the spin excitations is greatly reduced, 
and their bandwidth D„ ~ J becomes much smaller than 
Ep. A similar effect is known to occur in the strong inter- 
action limit of some lattice models, such as the Hubbard 
model. '^"'^ In this regime, the description of the spin exci- 
tations in the language of non-interacting spin waves is 



applicable only at very low temperatures T <C J. We 
show in this section that a generalized picture of decou- 
pled charge and spin excitations remains valid even at 
T > J. In this picture the charge excitations are still 
given by the waves of charge density (plasmons), and the 
spin waves are replaced with the excitations of a Heisen- 
berg spin chain. 



A. Luttinger-liquid picture of one-dimensional 
electrons systems 

The problems involving low-energy properties of inter- 
acting ID electron systems are conveniently described in 
the framework of the bosonization technique. ^^'^^ The 
first step in this approach is to linearize the spectrum 
of electrons near the Fermi level, thereby replacing the 
quadratic dispersion law e(fc) = ?i^fc^/2m with the lin- 
ear one. In this Tomonaga- Luttinger model the electrons 
are separated in two branches, the left- and right-movers, 
with energies eL,R{k) = hvpiTk — kp), where vp and kp 
are the Fermi velocity and Fermi wavevector. One can 
then present the fermionic field operators V'l.a and ipR,\ 
in terms of fields 0a and 9x satisfying bosonic commu- 
tation relations [(I)\{x),dy9\'{y)] = iT:5(x — y)6\x using 
the following rule 

■<PlAx) = ^4^e-^'='^^e*'^^(^)~*^^(^\ (la) 



V27ra 



(lb) 



Here A =tii is the spin index, a is the short distance 
cutoff, r\L.\ and r\i^^\ are Majorana fermion operators. '^^ 
In terms of the bosonic variables the Hamiltonian of 
an interacting ID electron system takes the forr-'^^^ 



(2) 



where the two terms Hp and describe the excitations 
of the charge and spin degrees of freedom, respectively, 
and have the forms 

H, = J^[n'KpIll + K;\d.<l>pr]dx, (3) 



2gi^ 
(27ra)2 



dx. 



(4) 



Here the new bosonic fields (pp^^r — {4'^ i (/)|)/\/2 and 
lip, cr = dxid] ± 6i)/tt\/2 satisfy the standard commuta- 
tion relations [0c((a;), Hq' (y)] — iS{x — y)Saa' and rep- 
resent the excitations of the charge and spin degrees of 
freedom. The Hamiltonian @-(@J) depends on five pa- 
rameters determined by the interactions between elec- 
trons: velocities Up and Ua- of the charge and spin ex- 
citations, dimensionless parameters Kp and K„, and the 
matrix element gi± of spin-flip scattering of a left-moving 
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electron and a right-moving one. In the absence of inter- 
actions, Up = Ucr = vp, Kp = Kg. = 1, and gi± = 0. The 
bosonized Hamiltonian correctly describes the charge and 
spin excitations of the system with energies below the re- 
spective bandwidths Dp^ ^ finup^a. 

In the most interesting case of repulsive interactions 
Kp < 1. The coupling constant gi±^ is positive and scales 
to zero at low energy scales D, 



gi± 



911- 



gi± 



In^ 
111 ^ 



(5) 



The parameter renormalizes along with 51^, ap- 
proaching the value Kg — 1 required by the SU(2) sym- 
metry as Ka = 1-1- 5ij_/27rMo-, Ref. Is^IbS 

Since the sine-Gordon term in Eq. Q vanishes at 
D/Da 0, the Hamiltonian ||2J) becomes quadratic and 
describes a Luttinger liquid.^- The latter represents a sta- 
ble fixed point of the problem, so the description of the 
system based upon the Hamiltonian ©-Q is expected 
to be valid in a broad range of interaction strengths. It 
is not immediately obvious, however, that the above pic- 
ture is applicable at low electron density n, i.e., when 
the interactions are so strong that the electrons form a 
Wigner crystal. Indeed, a Fourier expansion of the wave- 
function of an electron localized in a small region of size 
a <^ near a given lattice site involves the wavevectors 
in a broad range Sk ^ 1/a ^ n ^ kp- Thus the stan- 
dard procedure of linearization of the electronic spectrum 
near kp leading to the Tomonaga-Luttinger model is not 
justified in this case. In addition, each electron is con- 
structed out of waves with both positive and negative 
wavevectors, and the picture of two separate branches of 
left- and right-moving particles is not applicable to the 
Wigner crystal. We show now that even though the con- 
ventional derivation leading to Eqs. ©-I® is not justi- 
fied, this Hamiltonian does describe the low-energy prop- 
erties of a ID Wigner crystal. 



B. Charge and spin excitations in a Wigner crystal 

At low electron density, nuB ^ 1, the properties of 
the system are dominated by the Coulomb repulsion, and 
the electrons occupy fixed positions on the Wigner lat- 
tice. The first correction to this picture is due to the 
small vibrations of the lattice, analogous to phonons in 
conventional crystals. In the long-wavelength limit these 
phonons can be described in the framework of elasticity 
theory. In this approach the crystal is viewed as an elastic 
medium. Its motion is described in terms of the displace- 
ment u(x) of the medium at point x from its equilibrium 
position and the momentum density p{x) . The energy of 
the system can then be written as a sum of kinetic and 
potential energies, 



H = 



2ran 



dx. 



(6) 



The second term here is ^{d'^E/dn'^){Sn)'^, where E is 
the energy of the resting medium per unit length, and 
the density perturbation Sn is proportional to the defor- 
mation of the medium, Sn = —n d^u. The parameter 
s — {n/m){d'^E/d-n?) has the meaning of the speed of 
density waves (plasmons) in the Wigner crystal. 

The speed of plasmons in a ID system with true 
Coulomb interactions between electrons diverges in the 
limit of long wavelength. In practice, however, the inter- 
actions between electrons are usually screened at large 
distances by a remote metal gate. In the model where 
the gate is a conducting plane at a large distance d ^ n^^ 
from the Wigner crystal, the speed of plasmons is 



em 



In(Cnd), 



(7) 



where 8.0, Ref.lM 

The classical Hamilton function © can be quantized 
by imposing commutation relations [u{x) , p{y)\ = ihS{x— 
y). The resulting Hamiltonian describes the propaga- 
tion of the electron density excitations in a Wigner crys- 
tal, and is completely analogous to the term Hp in the 
Hamiltonian ©"(EJ of the Luttinger liquid. Comparing 
the commutation relations of the bosonic fields entering 
Hamiltonians and JBJ , and taking into account the ex- 



pressions for the density perturbation dn 
—ndxU, we identify the fields as 



V2 



u{x) = — <t)p{x), p{x) = ^ 
Tin y z 



Hp(x). 



(8) 



Using these expressions we can relate the parameters in 
the Hamiltonian Hp to the properties of the Wigner crys- 
tal as follows 



K 



(9) 



Here vp = iThn/2m is the Fermi velocity in a non- 
interacting Fermi gas of density n. 

The electrons in a Wigner crystal are repelled from 
each other by strong Coulomb forces. In the harmonic 
chain approximation we used so far the electrons are al- 
lowed to move about their equilibrium positions; how- 
ever, the amplitude of these oscillations remains small. 
As a result the electrons never move from one site on the 
Wigner lattice to another, and can be viewed as distin- 
guishable particles. Therefore the energy of the Wigner 
crystal state in this approximation does not depend on 
the electron spins. 

To account for the spin dependence, one has to in- 
clude the processes in which electrons tunnel through the 
Coulomb potential repelling them. Considering a pair of 
electrons at two neighboring sites of the Wigner lattice, 
one notices that depending on their total spin, the two 
electrons occupy either a symmetric or an antisymmet- 
ric state in the respective double-well potential. Thus 
the energy of the pair contains a term J Si ■ S2, where 
J > is the difference of energies of the antisymmetric 
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and symmetric states, and Si are the operators of elec- 
tron spins 1^ Taking into account all the nearest neigh- 
bor sites, we find that the spin properties of the Wigner 
crystal are described by the Hamiltonian of an antiferro- 
magnetic Heisenberg spin chain 



i+i- 



(10) 



Since the exchange is due to tunneling, the constant J is 
exponentially small, 



J — J* exp — 



V 



y/noB 



(11) 



To accurately evaluate J one has to take into account the 
fact that when two neighboring electrons tunnel through 
the Coulomb barrier separating them, all other electrons 
also move. Hausler"^^ suggested an approximation that 
neglects the motion of other electrons; his result corre- 
sponds to the value of 77 « 2.87 in a finite chain of 15 
electrons. We solve this model in the case of an infinite 
chain in Appendix 1X1 and obtain the value of 77 « 2.82. 
We also estimate the prefactor as 



J* w 1.79 



(nas)3/4' 



(12) 



where Ep = (7r?in)^/8m is the Fermi energy of a non- 
interacting electron gas of density n. 

The spin part (|10|l of the Hamiltonian of a ID Wigner 
crystal is very different from that of a weakly interacting 
electron gas, Eq. Q). It is easy to show, however, that at 
low energies D <ti J the two Hamiltonians are equivalent. 
To accomplish that we use the standard procedura^^''^^ of 
bosonization of spin chains. The first step is to perform 
the Jordan- Wigner transformation 



ajai~-, S'f +iSf = aj exp 



i-i 

•E 



a-ja-j I , (13) 



which expresses the spin operators in terms of creation 
and destruction operators and a of spinless fermions. 
In terms of these operators the Hamiltonian IjlOl) becomes 

Ha = ^'^J [ajai+i + aj^-^^ai^ 

+2 [ajai + + ■ (14) 

Thus the Heisenberg spin chain (|10|l is equivalent to the 
model (|14() of interacting lattice fermions. 

The second step is to bosonize the Hamiltonian (|14|l . 
At low energies one can replace the lattice model (|14f) 
with a continuous one, a/ — > a(jj), linearize the spectrum 
of the fermions near the Fermi level, and then apply a 
bosonization transformation 



0'L,R{y) 



(15) 



The resulting bosonized Hamiltonian of the spin chain is 
equivalents^ to Eq. I^J . The value of the speed Ua of the 
spin excitations is easily deduced from the Bethe ansatz 
solutionSL^ of the Heisenberg model, 



ttJ 



(16) 



Thus we have established that the bosonized Hamilto- 
nian lO^I© adequately describes the low-energy proper- 
ties of not only weakly interacting electron systems, but 
also of a ID Wigner crystal state at nas 'C 1. However, 
it is important to keep in mind that the applicability of 
the bosonized description to the Wigner crystal is limited 
to very low temperatures T <^ J. Given the exponential 
dependence Hll|l of the exchange constant on density, this 
condition can be easily violated even at fairly low temper- 
atures. In this case one has to use the more complicated 
form H10|) of the Hamiltonian H^. We show in Sec. El 
that this breakdown of the Luttinger liquid picture gives 
rise to significant deviations of conductance of quantum 
wires from the quantized value 2e^//i. 



C. Spin-charge separation at ultralow electron 
densities 

The Wigner crystal picture discussed in Sec. Ill Hl relies 
on the long-range nature of the Coulomb interaction po- 
tential V{x) — e'^/e\x\. In general, a ID electron system 
forms a Wigner crystal state at n — > only if the inter- 
action potential decays slower than l/jxp at x — > 00. In- 
deed, for potential V{x) oc l/|x|^ the interaction energy 
of two electrons at the typical interparticle distance 
is V oc rf , whereas the kinetic energy Ep oc n^. Thus at 
7 > 2 the interaction energy is negligible at n 0. 

The electron density in quantum wires is usually con- 
trolled by applying voltage to metal gates. The pres- 
ence of a gate affects the electron-electron interactions 
at large distances. For instance, if the gate is modeled 
by a conducting plane at a distance d from the wire, the 
interactions between electrons become 



{2df 



(17) 



The screening of the Coulomb potential by the gate re- 
duces the potential ((T7|l to V{x) — 2e^d^ le\x\^ at large 
\x\. Therefore, the Wigner crystal picture fails in the 
limit n — !■ 0. Comparing the interaction potential at 
the interparticle distance V{n^^^ with the Fermi energy 
of electrons, one concludes that within the model H17() 
the Wigner crystal state exists only in the density range 

As long as the system is in the Wigner crystal state, its 
spin excitations are described by the Heisenberg model 
(|10|l . However, the expression Hll|) for the coupling con- 
stant J relies on the pure Coulomb interaction between 
electrons. In the case of interaction potential screened by 
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the gate the exponential decrease of J with decreasing 
density stops at n ~ d~^, because the potential screened 
by the gates falls off rapidly at distances a; 3> d. Using 
the method described in Appendix IXI one estimates 

J^Ef[ exp -ryJ— (18) 

\aB J y S as j 

a,t as / (P ^ n ^ d^^ . In the case of interaction potential 
((T7|l the constant fj « 8.49. 

The distance to the gate in quantum wire devices is 
typically large, d ^ IOob, and most experiments are 
performed at densities well above as/d^. However, if 
the density is reduced to n < asld^, the Wi gner crys- 
tal picture used in Sec. IIIBI will fail. It is interesting 
to explore to what extent the conclusions of Sec. Ill Bl 
will be affected. To this end, let us now study the limit 
nd'^/as — > 0. 

At n <C as/d^ the interaction between two particles 
at a typical distance is small compared to their 
kinetic energy ~ Ep- On the other hand, when the 
distance between electrons is sufHciently short, \x\ ^ 
rr^{nd'^ / qbY^^ ^ ^ they experience strong repulsion 
V{x) Ep. Thus in the limit of low electron density one 
can model the interaction potential (|17|l by short-range 
repulsion 

V{x) = V5{x). (19) 

The constant V should be chosen in such a way that the 
scattering phase shift for two electrons interacting with 
potentials V{x) and V5{x) are identical. For the model 
((T7|l this condition gives 

The exponentially large value of V reflects the fact 
that the strong repulsion H17|) leads to almost perfect 
backscattering of electrons off each other. 

At V ^ oo the electrons are separated by thin hard- 
core potentials. In this limit they can be viewed as dis- 
tinguishable particles, and the eigenvalues of energy be- 
come independent of the electron spins. The wavefunc- 
tions of the system essentially coincide with the Slater 
determinants for spinless non-interacting fermions. Upon 
bosonization, the Hamiltonian Hp of this system takes 
the form ©■ The plasmon velocity s in this system is 
the Fermi velocity of non-interacting electron gas of den- 
sity n, which is twice the Fermi velocity of non-polarized 
electron gas, s = 2vf- Thus according to Eq. (j5J we 
haveSi Kp = 1/2. Additional properties of this model 
were recently discussed in Refs. 39,40. 

At large finite V the electrons can change places as a 
result of scattering, and the energy acquires a weak de- 
pendence on the spins. This dependence can be deduced 
from the well-known properties of the one-dimensional 
Hubbard model. It has been shown by Ogata and 



Shiba''^''*^ that at U /t oo the spin and charge exci- 
tations of the Hubbard model are completely separated, 
with the Hamiltonian of spin excitations taking the form 
of the Heisenberg model 110|) . The magnitude of the ex- 
change constant in this Hamiltonian was found^ to be 

4<^ / sin27rne\ / n 

where is the average number of electrons per site. In 
the limit rig — > the Hubbard model is equivalent to an 
electron gas with quadratic spectrum and point-like in- 
teraction (|19|l . The limiting procedure can be performed 
by introducing infinitesimal lattice period a in the Hub- 
bard model, identifying the parameters t = ?i^/2ma^, 
U — V /a, Tie — na, and taking the limit a — > 0. Apply- 
ing this procedure to the formula H21() . we find 



Using the estimate ()20|) of parameter V for the interaction 
potential (fTTji . we find 

j^Ep—e^^{~fl.p\. (23) 
as \ \ as j 

Note that the our results H18|l and H23|l for the ex- 
change constant are of the same order of magnitude at 
n = as I iff' ■ 

So far we have demonstrated that the description of 
the system in terms of the Hamiltonian in spin-charge 
separated form H = Hp + H„, with Hp and H„ given by 
Eqs. and (|10() is valid in two different regimes. The 
first one is the Wigner crystal state at electron densities 
in the range as /d"^ <^ n a~^^ , and the second is the low 
density limit n <C as/d'^, where the picture of point-like 
interactions (|19|l is applicable. One can show*'^ that in 
fact this picture of spin-charge separation holds at any 
density n <^ a^^, including the regime n ~ as/d^. 

The exchange constant J in the effective spin chain 
Hamiltonian H10|l monotonically decreases as the electron 
density n is lowered. In the most interesting range of 
densities d~^ <C n <C a^^ the dependence of exchange 
on n is exponential, Eq. Hll() . At lower densities the 
dependence becomes a power-law one. Specifically, in 
the density ranges as/d^ ^ n <^ d^^ and n <C as/d^ 
one can use the estimates (|18|l and (|23|) . respectively. 

For the sake of simplicity, in the following sections 
we assume that the electron density is in the range 
as/d^ ^ n <C a^^, and refer to the electron system 
as a Wigner crystal. However, all of our conclusions re- 
main valid at any densities n <$i a~^^ , if the value of the 
exchange constant J is adjusted as discussed in this sec- 
tion. 
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III. CONDUCTANCE OF A QUANTUM WIRE 
WITH SPIN-CHARGE SEPARATION 



The spin-charge separation has a profound effect on 
the conductance of quantum wires. Indeed, the electric 
field applied to the wire couples to the electron charges 
and has no effect on spins. As a result, the spin degrees 
of freedom remain decoupled from charge ones, and the 
rather complex form of the Hamiltonian iJ^- has no effect 
on the conductance. In this section we review the known 
results for the conductance of a quantum wire with spin- 
charge separation. 



A. Infinite wire 

Conductance of an infinite Luttinger liquid is given by 
G — 2Kpe^/h. This result was obtaincd"''* '^^ by assuming 
that a weak electric field is applied to a small part of the 
wire, and the conductance was evaluated by using the 
Kubo formula. In the following sections it will be more 
convenient to evaluate the conductance of the Wigner 
crystal in the regime of applied current. It is therefore 
instructive to reproduce the result G = 2Kpe^/h in this 
approach. 

Let us consider a quantum wire whose charge dynam- 
ics are described by the Hamiltonian , and enforce the 
current / = Iq cos ujt at x = 0. By doing so we impose 
a boundary condition upon the charge field 0p(O,i). In- 
deed, the bosonization expression for the electric current 
is / = e-^(/>p. (In the case of a Wigner crystal, this can 
be checked by using Eq. JHJ and the definition / = enii 
of current in terms of the velocity ii of the crystal.) Thus 
the field chn satisfies the condition 



where the function 



■ sinLot. 



(24) 



(25) 



is related to the current as I — eq and has the meaning of 
the number of electrons that passed through point x = 
at time t. 

By imposing a time-dependent boundary condition 
(|24|l we drive the system with an external oscillating 
force. This leads to emission of plasmon waves and dis- 
sipation of the energy from the driving force to the infi- 
nite Luttinger liquid. We will find the resistance of the 
wire Rp by evaluating the energy W dissipated in unit 
time and comparing the result with the Joule heat law 
W = hl^Rp. We present a formal derivation in Ap- 
pendix IB] here we limit ourselves to a simple semiclassi- 
cal argument. 

Solving the Hamilton equations with Hamiltonian ^ 



and boundary condition H24|l we findS 
tt/o 



(t)p{x,t) 

np[x,t) 



\/2euj 
lo 



V2eKp 



smtoit — \x\/up), (26a) 
cosuj{t — \x\/up). (26b) 



Substituting this solution back into Eq. (|3J), we find the 
following expression for the time-averaged energy density 
in the Luttinger liquid. 



(27) 



The plasmon wave (1^ carries the energy {E)t at speed 
Up in two directions. Thus the total energy dissipated 
into plasmon waves in unit time is given hyW — 2up{E)t. 
Comparing this result with W = i/p-Rp, we find the 
resistance 



Rn 



2Kpe'^ ' 



(28) 



in agre ement with the result for the conductance found 
in Rcf . I44l45i 



B. Finite-length quantum wire between two 
non-interacting leads 

The result G = 2Kpe^/h indicates that in a quan- 
tum wire with repulsive interactions conductance should 
be below the quantized value 2e^/h. Furthermore, it is 
expected to decrease as the electron density n is low- 
ered. However, the experiments consistently show per- 
fect quantization^^ of conductance at 2e'^/h in a broad 
range of n. 

This controversy was resolved^Siiiii£ by noticing that 
instead of an infinite quantum wire, the experiments 
study transport through a finite-length wire connecting 
two bulk leads. Since the leads are not one-dimensional, 
their properties are not adequately described by the Lut- 
tinger liquid model @-l0J. Instead, the electrons in the 
leads are expected to be in a Fermi liquid state. 

To find the conductance of such devices, one can 
j-^Q jgj46j47j48 ^j^g leads connected to the wire by two semi- 
infinite non- interacting wires. In this model the system 
remains one-dimensional, but the interactions are non- 
vanishing only in the central part of the system. The 
length L of the interacting part is identified with the 
length of the wire. Assuming that the interactions fall 
off gradually at a; — > ±oo, one can neglect the back- 
scattering of electrons from the interacting region. In 
this limit the charge dynamics are still described by the 
Hamiltonian ||2J), but the parameters Kp and Up become 
functions of coordinate x. 

The measurements of dc conductance in experiments 
are conducted at very low frequencies lu <JC Up/ L. Thus 
the wavelength of the plasmons emitted in the system is 



7 



much greater than the length of the interacting region 
L. Consequently, one should use the value of the pa- 
rameter Kp in Eq. H28|l taken at x ^ ±00, i.e., in the 
non-interacting region, where Kp = 1. Thus the resis- 
tance of the device becomes 



restoring the perfect quantization. Careful treat- 
ments^Si^ii^ of the problem lead to the same conclusion. 

IV. VIOLATION OF SPIN-CHARGE 
SEPARATION IN QUANTUM WIRE DEVICES 

As we saw in Sec. lIIII the inhomogeneity of the system 
caused by coupling of the wire to the leads changes the 
conductance from 2Kpe^/h to 2e'^/h. This conclusion 
was derived from consideration of the charge excitations 
only, as the spin degrees of freedom were assumed to be 
completely separated. We now turn to the effect of the 
inhomogeneity on the spin part iJ^. of the Hamiltonian. 

We assume that the central part of the wire contains 
a purely one-dimensional electron system at low density 
n <C ttg^, so that the Wigner crystal model is appropri- 
ate. The wire is also assumed to be smoothly connected 
to the leads, where the effective interactions are weak. 
This is due to several effects. First, the electron density 
grows as one moves away from the wire into the leads. 
This effectively reduces the interaction strength, as the 
latter is characterized by parameter (nos)"^. In addi- 
tion, the wire becomes wider when it couples to the leads. 
As a result, when two electrons arrive at the same coor- 
dinate X along the wire, they are no longer as close to 
each other as in the middle of the wire. This reduces the 
strength of interactions between ID electrons. The two 
mechanisms have very similar effect on the Hamiltonian 
H^. For simplicity, in the following we limit our discus- 
sion to the effect of inhomog eneous el ectron density. 

Following the ideas of Refs. we model the wire 

connected to the leads by an inhomogeneous ID system. 
The main source of inhomogeneity is the dependence n{x) 
of the electron density on position. We assume that the 
density takes a constant value n{x) — n inside the wire, 
i.e., at |a:| < L/2, and gradually grows to a very large 
value rioo ^ at a; ^ ±00. 

In experimental devices the dependence of electron 
density on the coordinate along the wire is caused by 
inhomogeneity of the external confining potential. Apart 
from changing the electron density, the external potential 
may also lead to backscattering of electrons in the wire. 
In a sufficiently long wire such processes may greatly sup- 
press the conductance at low temperature 4^*^ In the 
Wigner crystal picture this phenomenon is interpreted 
as pinning of the crystal by the external potentiali^ 
On the other hand, the best available experiments show 
good quantization of conductance, indicating that the 
backscattering remains negligible. This is most likely the 



result of smoothness of the confining potential. Indeed, 
the backscattering involves the change of the electron 
wavevector by 2kF- Thus an external potential that is 
smooth at the scale of interparticle distance will 
cause exponentially weak backscattering. In this pa- 
per we assume that the external potential is sufficiently 
smooth, so that the backscattering can be neglected. 

Under the above conditions the low-energy properties 
of the system may be described by the bosonized Hamil- 
tonian 0-101), but with position-dependent parameters 
Mp,£r, Kp^a-, gi±- In this paper we assume that the temper- 
ature is small compared with the bandwidth Dp ^ hnup 
of the Hamiltonian Hp, so that the discussion of the ef- 
fect of the charge modes on the conductance presented in 
Sec. lIIIl is valid. On the other hand, we will be interested 
in the case of temperature comparable with the band- 
width Z?CT J of the Hamiltonian H^. In this regime the 
bosonized version of He is not applicable, and one 
should instead use the Heisenberg model H1U|I . 

Since the exchange constant 1)11(1 strongly depends on 
the electron density n{x), the parameter J in Eq. H1U|) 
should also be considered position-dependent. In partic- 
ular, the strength of the exchange coupling between the 
two spins at the neighboring sites I and Z -1-1 of the Wigner 
lattice is a function of the coordinate xi of the l-th. elec- 
tron: J — J{xi). It is important to note that in the 
presence of electric current / the Wigner lattice moves, 
so the coordinate xi of the l-th lattice site depends not 
only on /, but also on time. 

The time dependence of xi can be accounted for by 
noting that if during the time interval t a number q{t) of 
electrons have moved from the left lead to the right one, 
the Z-th site of the lattice has shifted to the {I + q)-th. 
position. Thus the time dependence of the positions of 
the lattice sites can be accounted for by replacing I 
I + q{t), and the Hamiltonian TJ^ takes the form 

He=J2j[l + q{t)]Si- Si+i. (30) 

Note that in this approximation the electric current / = 
eq(t) is assumed to be uniform throughout the wire. This 
is true in the dc limit uj <C Up/ L. 

It is important to note that the form (|30|l of the Hamil- 
tonian Ha- violates the spin-charge separation. Indeed, 
the coupling between the spins depends on the amount 
of charge that passed through the wire, which is related 
to the field 0p, see Eq. (1^ . As a result, the conduc- 
tance of a quantum wire connected to bulk leads may be 
affected by the spin excitations.^'' 

To find the effect of spin subsystem on the conduc- 
tance, one could substitute the expression (|24(l for q{t) 
into Eq. (|30|l . and consider the complete Hamiltonian 
Hp + Ha- without relying on spin-charge separation. In 
this approach one needs to add to the Hamiltonian a 
term describing the applied bias, and evaluate the elec- 
tric current. However, it is more convenient to treat the 
current I(t) in the wire as an external parameter. In this 
case q(t) is also a parameter, and the Hamiltonians Hp 
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and Her still commute. The only consequence of the vi- 
olation of spin-charge separation in this approach is the 
dependence of Her on the current I{t). 

The presence of an oscillating parameter (|25|l in the 
Hamiltonian (|30|l may lead to creation of spin excita- 
tions. Using the approach of Sec. IIIII we will calculate 
the energy dissipated into spin excitations in unit time. 
In the limit of weak current, the dissipation is found in 
the second order of the perturbation theory in the am- 
plitude Iq of current oscillations. Thus in addition to the 
plasmon result for the energy W dissipated in unit time, 
we will obtain a similar contribution of the spin modes: 



(31) 



Comparing this result with the Joule heat law W — 
^IqR, we conclude that the resistance R of the wire is 
given by sum of two independent contributions. 



R — Rp + Ra 



(32) 



Note that the first term in this expression is already 
known, Eq. (|29|l . The second term is discussed in Scc.lVl 
It is interesting to point out that the result (|32|l may 
be interpreted as a total resistance of the charge and 
spin subsystems connected in series, whereas naively one 
might expect a parallel connection. The reason is that 
the spins do not directly respond to the applied voltage, 
as required for the latter interpretation. Instead, the 
spin subsystem responds to the electric current. Thus 
the Hamiltonians Hp and H„ become independent in the 
regime of applied current, in analogy with the problem 
of two independent resistors connected in series. A result 
similar to Eq H32fl has been obtained for the resistivit y of 
two-dimensional strongly interacting systems in Ref. 




nL/2 nL/2 



FIG. 1: Scketch of the dependence J[y] in our model. Inside 
the wire, |y| < nL/2, the exchange J is exponentially small, 
Eq. IIH . As one moves toward the leads, J grows, and at 
J/ cx) it saturates at Joo ~ Ep. 



the exchange is exponentially small, Eq. As the 

wire connects to the bulk leads, the density n(x) begins 
to grow. At nas ^ 1 the exchange J becomes of the 
order of the Fermi energy, see Eqs. Hll() . 1)121) . 

Strictly speaking, the Wigner crystal picture is valid 
only at naB ^ 1, that is as long as J ^ Ep- On the other 
hand, we will be interested in the properties of the system 
at low energies D <C Ep. Thus at J ~ Ep when the 
Wigner crystal picture fails, we are only concerned with 
the energy scales much lower than J. As we saw in SecUll 
at those scales one can use the bosonized Hamiltonian 
Q regardless of the applicability of the Wigner crystal 
description. Thus we can ignore the difference between 
the Wigner crystal and weakly interacting electron gas 
at large density n ^ as, and simply assume that in the 
leads the exchange J saturates at J^o ^ Ep. 

The properties of the function J[y] can thus be sum- 
marized as follows. 



^^^^ (j-tiEp, at|?/|<nL/2, 
1 Joo Ep, at \y\ ±oo, 



(34) 



V. SPIN CONTRIBUTION TO THE 
RESISTANCE 

To find the contribution R^ of the spin subsystem to 
the resistance of the device, we study the dissipation of 
energy into spin excitations caused by the time depen- 
dence of the Hamiltonian H3U|I . We start by performing 
the Jordan- Wigner transformation (|13|l and converting 
the Hamiltonian to the fermionic form 

Ha = + lit)] (ajai+i + al_^^-^ai 



+2 (ajai + i j (aj^^ai+i + ^ 



(33) 



In the absence of the external magnetic field the average 
z-component of the spin at every site of the lattice must 
vanish. Thus according to Eq. (|13|l the occupation of 
each site is (ajai) = ^. This means that the Fermi level 
is in the middle of the band, /i = 0. 

The exchange J[y] strongly depends on the position y. 
Inside the wire the electron density is low, nas ^ 1, and 



see Fig. ^ Note that y is the coordinate on the Wigner 
lattice. Since we consider the limit of very smooth con- 
fining potential, all the physical quantities change very 
little at the interparticle distance. We therefore assume 
that J[y] is a slowly varying function: \dJ/dy\ <^ J[y]- 



A. XY model 

The Hamiltonian H33|) describes a system of strongly 
interacting fermions. As a first approximation we will 
simplify the problem by neglecting the interactions be- 
tween the fermions. 



(35) 



This Hamiltonian corresponds to the fermionized version 
of the XY model of a spin chain, in which the coupling 
of the z-components of spin operators is neglected. This 
approximation violates the SU(2) symmetry of the prob- 
lem, and is therefore rather crude. On the other hand. 
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the resistance Ra can be found exactly for model (|35|l . 
and the result will provide considerable insight into the 
properties of model (jSSll- 

Hamiltonian (|35|l represents an inhomogeneous version 
of the the tight-binding model of lattice fermions. In the 
uniform case, J[y\ ~ const, the spectrum is well known, 

e{k) = Jsinfc, (36) 

where the wavevector k is measured from kp — 7r/2. One 
can either assume that k varies in the interval — tt < fc < 
TT, or choose < A: < tt and treat Eq. (|36|l as spectra of 
two branches of excitations, the particles and holes. 

In the absence of electric current in the wire one can 
omit q[t) in the Hamiltonian 135|) and view it as a tight- 
binding model with slowly varying bandwidth 2J[?/]. In 
the leads the bandwidth 2Jao is very large; it narrows 
down to a very small value 2 J in the wire, Eq. (|34|l . 
The particles moving toward the wire in one of the leads 
cross to the other lead if their energies are below the 
small exchange J in the wire; the particles with e > J 
are reflected. 

In the presence of the electric current /, the constric- 
tion of the band in the Hamiltonian H35|l moves with re- 
spect to the lattice with velocity q — I /e. The particles 
reflected from the moving constriction change their en- 
ergy. These processes lead to the dissipation of energy 
and contribute to the resistance R^- 

For non-interacting fermions, the problem of evaluat- 
ing the energy W dissipated in unit time by a moving 
scatterer can be solved for arbitrary reflection coefficient 
7?,(e), see Appendix lO Here we find W in the semiclassi- 
cal limit, which is valid for very slowly varying bandwidth 
J[y], when 7^(e) = 9{e~ J). 

In the limit of slowly varying J[y\ one can apply the 
result for the spectrum of particles at every point 
in space, and treat the excitations as classical particles 
with energy 

H{y,p,t) = J[y + vt]sm^. (37) 

Here y is the coordinate of the particle, p is its momen- 
tum, and V = I/e is the velocity of the constriction. 
For simplicity we will consider the case of dc current, 
/ — const. To find the linear conductance of the quan- 
tum wire, one can limit oneself to the case of very small 
current, and assume v <^T /Ti, J /Ti. 

The time-dependent energy H37|l should be treated as a 
Hamilton function, and the trajectory of the particle can, 
in principle, be found by solving the classical Hamilton 
equations. One can easily check that the quantity 

E{y,p,t)^H{y,p,t)+pv (38) 

is an integral of motion. It has the meaning of energy 
of the particle in the frame moving at the speed of the 
constriction. 

A particle with energy e ~ T moving in the right di- 
rection has a very low momentum p when it is in the 



leads, p/fi = e/Joc ^ T/Ep ^ 1. Thus its integral of 
motion E{y,p,t) = e. As the particle approaches the 
constriction, its momentum increases, so that E retains 
its value despite the decrease of the bandwidth J. At 
small V the maximum allowed value of E in the wire is 
reached at p — TTh/2 and equals i?max = J + ^hv. Thus 
at e < J + ^hv the right-moving particle moves from 
the left lead to the right one, and its energy e remains 
unchanged. If the energy e exceeds J -I- ^hv, the parti- 
cle cannot enter the wire. When its momentum reaches 
7Tti/2 at a point to the left of the constriction, the particle 
is reflected. Deep in the left lead its momentum is very 
close to Trh. Due to the pv term in integral of motion 
its energy H decreases to e — irfiv. 

Similarly, since a left-moving particle in the right lead 
with energy e has momentum very close to nh, its integral 
of motion H38() is E = e + nhv. The condition E < i?max 
for transmission through the constriction for such parti- 
cles means e < J — ^hv. As the particle reaches the left 
lead, the momentum is again near ttH, i.e., conservation 
of E results in conservation of energy H — e. On the 
other hand, particles with energies e > J ~ ^hv are re- 
flected back to the right lead, and their momentum on 
the right-moving branch is near p = 0. Thus the energy 
of these particles increases from e to e -t- irhv. 

To summarize, the particles in the leads with energies 
e < J — ^hv cross the constriction region without change 
of energy. The particles with energies e > J + ^Tiv are 
always reflected by the constriction. The ones in the 
left lead decrease their energy by nhv, while the ones in 
the right lead increase their energy by the same amount, 
so that these contributions to the total energy of the 
system compensate each other. Finally, in the narrow 
range of energies J — ^hv < e < J + ^hv the right- 
movers go through the constriction without change of 
energy, whereas the left-movers reflect back to the right 
lead with energy gain nhv. 

The total current of left moving particles and holes in 
the narrow energy interval of width nhv near e = J is 
given by 5N = {2/h){nhv)f{J), where /(e) = l/(e<^/^ + 
1) is the Fermi function. Thus the total energy trans- 
ferred to the spin excitations in unit time is 

W'''' = nhv^fiJ)^l'^fiJ). (39) 

Comparing this result with the Joule heat law W = 
I^R, we obtain the spin contribution to the resistance 

RV = ^J{J)- (40) 

At low temperature T <ti J most of the particles have 
energies below J and pass through the constriction elas- 
tically. Only an exponentially small fraction of particles 
are reflected and contribute to the dissipation. Thus the 
result 1401) is exponentially small at low temperatures, 
jjjfy ~ (ft,/2e^)e~'^/^. As the temperature is increased, 
a greater fraction of the particles are reflected by the con- 
striction, and the R^^ increases. In the limit T/J — > oo 
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all the particles are reflected, and the resistance saturates 
at = h/Ae^. 

In this section we studied the simplified XY model, 
in which the z-component of coupling in the Hamilto- 
nian H^Ofl was neglected. Thus the result (|40|1 cannot be 
applied directly to the problem of conductance of a quan- 
tum wire in the Wigner crystal regime. However, much 
of the physics leading to Eq. (|40l) can be carried over to 
the case of the isotropic model (|5n|) . 



B. Isotropic coupling 

The problem of the isotropic inhomogeneous Heisen- 
berg spin chain (|30|l is far more complicated than that 
of the XY model (|35|l . However, it can still be some- 
what simplified by assuming that J[y\ is a very slowly 
varying function. Then each moderately long section of 
the spin chain can be approximated by the homogeneous 
Heisenberg model. The latter allows for exact solution^ 
by means of Bethe ansatz. The low-energy excitations 
of the isotropic Heisenberg spin chain are spinons with 
energy spectrur'^^^^ 



e(fc) — sin k. 



(41) 



Although the spinons do not obey Fermi statistics, the 
similarity between Eq. (|41l) and the spectrum H36|l of the 
excitations of XY model enables us to find the tempera- 
ture dependence of i?cr at T <C J. 

Indeed, most of the discussion leading to Eq. did 
not rely on Fermi statistics of the excitations. One can 
apply the arguments of Sec. IV Al to the problem of scat- 
tering of spinons by the constriction of the band in the 
wire. In particular, one concludes that spinons with en- 
ergies below TT J/2 pass through the constriction without 
scattering and do not change their energy. Thus the dis- 
sipation is exponentially small at T <C J, and one finds 
Ra (X exp(-7r J/2T). 

Since the occupation of states with high energy e = 
7rJ/2 at low temperature is exponentially small and in- 
dependent of statistics, on can naively expect the resis- 
tance Ra to be given by the low-temperature asymptotics 
of Eq. H40|l upon replacement J — > 7rJ/2. Then one ob- 
tains 



, ttJ 

R„ = i?oexp ( - — 



(42) 



This approach gives the prefactor Rq — h/2e^. 

Unfortunately the analogy between spinons and 
fcrmion excitations of the XY model does not enable one 
to find the prefactor in Eq. (|42|1 . Unlike the excitations 
of the XY model, the spinons interact with each other, 
and the energy of a spinon is affected by the presence of 
other spinons. In the limit T ^ the density of other 
spinons is small, and the energy is given by Eq. H41|) . At 
finite T the result 141|) may acquire a small correction. 



The exponent of Eq. H42() is determined by the maximum 
energy of a spinon tt J/2. Even a small correction to this 
energy may affect the prefactor Rq. 

At high temperature T ^ J the resistance contribu- 
tion R^ evaluated within the XY model approximation 
saturates, because in this regime all the excitations are 
refiected by the constriction. This feature is preserved in 
the model with isotropic coupling, as at J ^ the spin 
excitations cannot propagate through the wire. The sat- 
uration value of i?o- at high temperatures can be found 
by noticing that in the part of the system away from 
the constriction, where J[l] ~ J^o ^ T, one can still 
bosonize the Hamiltonian (|33() and use the form In- 
side the constriction the bosonization is not applicable, 
and this region is modeled by imposing a boundary con- 
dition on the bosonic field 0o- corresponding to the fact 
that there is no spin propagation between the regions of 
large positive and negative I. 

The form of the boundary condition can be deduced by 
replacing the constriction of the bandwidth ^ 
J[l] in the Hamiltonian H33|l with a high potential bar- 
rier for the fermions. The barrier is modeled by a large 
backscattering term u(a^afl + a^a^,) at site I ~ ~(l{t)- 
Upon the bosonization transformation (|15|l this term be- 
comes —V cos[\/2<j)cr{y) — 2kFy]\y=~q{t), with kp on the 
lattice being tt/2. Since this scattering term is very 
strong, V CO, it pins the field (l)a{—qit),t) to the 
value —y/2kFq{t) ~ — (7r/\/2)(z(t). This time-dependent 
boundary condition leads to the emission of spin waves, 
in analogy with Sec. lIII ^ where the boundary condition 
(|24|l gave rise to plasmon waves (I26f) . In the limit of weak 
current, / ~ eujq — *■ 0, the wavelength of the spin waves 
^ Joo/fiuj is much larger than q{t), and instead of impos- 
ing the boundary condition at y = —q{t) one can impose 
it at ?/ = 0. Then the boundary condition becomes 



'V2 



q{t). 



(43) 



Note that up to inessential negative sign Eq. H43|) is 
identical to the boundary condition (|24|l . The respective 
Hamiltonians JSjl and Q are also essentially identical at 
low energies, as the sine-Gordon term is irrelevant. One 
can therefore carry over the results of Sec. IIII Al for the 
dissipation of energy into plasmon waves and the result- 
ing contribution to the resistance. Adapting Eq. (|28|) to 
the parameters of Hamiltonian Q , we find 



R<j 



(44) 



In the dc limit the frequency of the driving force — > 0, 
and the wavelength of the spin waves is very long. Thus 
the parameter Kcr in Eq. (|44|l is taken at large distances 
from the constriction, where the SU(2) symmetry de- 
mands = 1. Consequently the spin contribution 
to the resistance in the model with isotropic coupling 
is given by 

h 

2^' 



Ra 



(45) 
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On the other hand, the XY model H35|l does not pos- 
sess the SU(2) symmetry, and the bosonization proce- 
dure (|15|l gives the quadratic part of Hamiltonian Q 
with K„ ^2. Then Eq. (O predicts R^^ = h/Ae^, in 
agreement with T J asymptotics of Eq. . 



VI. DISCUSSION OF THE RESULTS 

The quantity most commonly measured in experiments 
with quantum wire devices is the linear conductance. In 
our theory its value is given by 



G 



1 



Rn + Ra ' 



(46) 



c.f. Eq. The contributions Rp and R^ to the resis- 

tance of the wire are determined by the properties of the 
charge and spin excitations of the system, respectively. 

Throughout this paper we consider the case of rela- 
tively low temperature, T <^ Dp 
the contribution of the charge 
h/2e'^ (see also Sec. UlTBI) 



R, 

ture above Dp leads to thermal smearing of conductance 
plateaus. No interesting electron correlation effects are 
expected in this case. 



--^ hnup. In this regime 
modes is well known: 
Raising the tempera- 



D 



At not too low electron density n 
D„ 



> 



the bandwidth 



Ep is the only relevant energy scale of the 
problem. Then at T <C Ep the contribution R„ van- 
ishes, and the conductance takes the well-known quan- 
tized value G = 2e^//i. On the other hand, in the in- 
teresting case of low density n <C a^^ another energy 
scale, the exchange constant J, appears in the problem. 
This scale is exponentially small, Eq. (|ll|l : in particular, 
J ^ Dp. In the limit of low temperature T ^ the con- 
tribution still vanishes. More specifically, at T ^ J 
we predict activated temperature dependence H42|) of R^, 
with activation temperature 'kJ/2. At higher tempera- 
tures Ra grows, and at T ^ J it saturates at the uni- 
versal value = /i/2e^, see Eq. (|45f) . Combining these 
results with Eq. we obtain 



h 



(47) 



This is our main result. It corresponds to an additional 
quantized plateau of conductance of a quantum wire at 
low electron density. The value of the conductance at this 
plateau is exactly one half of the quantized conductance 

The plateaus of conductance at e'^/h have 
been observed at low electron densities in several 
experiments.ASJiiSJi The authors of Refs. Iinlllll2ll3l 

attributed this feature to the spontaneous spin polariza- 
tion in quantum wires. This interpretation contradicts 
to the theorem by Lieb and Mattis;^ stating that 
the ground state of a ID electron system cannot be 
spin-polarized in the absence of magnetic field. One 
can hypothesize that the ferromagnetism in quantum 



wires is possible because the electrons are not truly 
one-dimensional; however to the best of our knowledge, 
no such theory is available at this time. In our theory the 
spin structure of the Wigner crystal state is described 
by the Heisenberg model (|10|l with positive exchange 
constant J, corresponding to antiferromagnetic coupling. 
Thus the ground state of the Wigner crystal in not 
spin-polarized, in agreement with the theoremn^ 

The temperature dependence of the conductance of a 
quantum wire device obtained in this paper is similar to 
the behavior observed in experiments on 0.7 structure in 
quantum point contacts i^s^i^ In agreement with exper- 
iments, conductance (146(1 remains 2e^/h at T ^ 0, but 
develops a negative correction at finite temperature. The 
activated temperature dependence of the correction fol- 
lowing from Eq. (|42|l is consistent with the measurements 
of Ref. 0. At high temperature the correction saturates, 
and the conductance develops a new plateau. Contrary 
to the experimentSf^iSii this plateau is at one half of the 
quantized value 2e^/h, rather than at 0.7 x (2e^/ft.). The 
relation between the plateau at e^/h and the 0.7 struc- 
ture was studied experimentally in Ref. Il2l It was found 
that the quasi-plateau at 0.7x (2e^/ft.) is observed in short 
wires, whereas in longer wires it shifts toward e'^ /h. In 
this paper we assume that the wire is long, so that the 
parameters of the system, such as the confining potential, 
Fermi energy, and exchange constant J, do not change 
significantly at the interparticle distance. It would be in- 
teresting to generalize our approach to the case of shorter 
wires and see whether the physics discussed in this pa- 
per may be responsible for the 0.7 structure in quantum 
point contacts. 

To test the relevance of our theory to the 
experimentsiSiiiiiSii^ showing plateaus at e^/h, one can 
check whether the experimental temperature exceeds the 
exchange energy J. Due to the strong exponential depen- 
dence (|ll|l of J on the density, the uncertainly of n may 
make the estimation of J difficult. Instead one may be 
able to determine J experimentally by applying magnetic 
field. Indeed, if the magnetic field B exceeds a certain 
critical value Be oc J, the spin chain becomes completely 
spin polarized. The magnitude of the critical field Be 
can be found by considering the spin-polarized state in 
a strong field B with a single spin-flip excitation. The 
energy of such an excitation \g\fiBB is reduced by 2 J due 
to coupling to neighboring spins. (Here g is the Lande 
factor, and fj,B is Bohr magneton.) Thus the complete 
polarization occurs at B > Be, where 



Br 



2J 



(48) 



By measuring the critical field Be required to achieve 
complete polarization of the spin chain one can determine 
the exchange constant J. 

In the case of non-interacting electrons at zero temper- 
ature the conductance does not depend on the magnetic 
field and remains 2e^ jh until the electron gas becomes 



completely spin polarized a,t B > B, 



(0) 



Ep/MglfiB- 
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In a polarizing field only one spin channel is allowed 
in the wire, and the conductance reduces to e^jh. In 
the case of a quantum wire at low electron density this 
behavior is preserved, but the step in conductance oc- 
curs at a much lower critical field (|48|l . Indeed, although 
in the presence of magnetic field spinous are no longer 
the elementary excitations of a spin chain, at B < Be 
one can introduce modified elementary excitations with 
similar properties. Then by repeating the arguments 
of Sec. IV Bl one concludes that the low energy excita- 
tions present in the system at T ^ cross the wire 
elastically, resulting in Ra- = and total conductance 
G = 2e^//i. At B > Be the wire is completely spin po- 
larized, and the spin excitations in the leads are reflected 
by the wire. This situation is completely analogous to 
the case of high temperature considered in Sec. IV Bl In 
particular, the resistance can be found by bosonizing 
the electron system in the leads and imposing the bound- 
ary condition (|43|l on the field This again leads to 
— ft./2e^ and reduces the conductance of the device 
to e^/h. Thus one can find the critical field (|48|l and the 
exchange constant J by measuring the magnetic field at 
which the conductance drops from 2e^//i to e^/h. 

Apart from the experiments with GaAs quantum wires, 
quantization of conductance at G — e^/h in the absence 
of magnetic has been observed in carbon nanotubes.-^ 
This anomalous quantization occurs when the current is 
forced to flow through the narrow tip of the tube. At 
small radius of the nanotube the Coulomb interactions 
between electrons become effectively stronger, and could 
conceivably suppress the exchange coupling J of the elec- 
tron spins below the temperature. Our result (|47|l would 
then explain the experimental data^ 
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APPENDIX A: ESTIMATE OF THE EXCHANGE 
CONSTANT J 

Here we estimate the exchange constant in an infinite 
ID Wigner crystal with the lattice constant b = 1/n. Fol- 
lowing the idea of Hausler— we evaluate J for two spins 
at neighboring sites I = and I = 1 using an approxi- 
mation where the only dynamical variable is the distance 
X ^ xi ~ xo between the two electrons. In this approxi- 
mation, xq + xi = b and all the other electrons (/ ^ 0, 1) 
are at fixed positions xi = lb. Then the Coulomb poten- 



tial takes the form 



Uix) 



e\x 




(Al) 



This potential has two degenerate minima at x = ±b 
corresponding to xq — 0, Xi — b and Xq — b, xi — 0. 
In the limit of strong Coulomb potential the tunneling 
between these two states gives rise to exponentially small 
splitting AE of the doublet. 

The ground state wavefunction of this system is an 
even function of x, and is therefore symmetric with re- 
spect to permutation xq <-> xi, while the first excited 
state is antisymmetric. The two states correspond to the 
values of the total spin of the two-electron system S = 
and S = 1, respectively. Thus the energy of the two 
components of the doublet can be written in terms of the 
electron spin operators at the two sites as Eq + J Sq ■ Si, 
where J is identified with the level splitting AE. 

Strictly speaking the infinite series in Eq. (HU di- 
verges. This is due the long range nature of the Coulomb 
interactions. In practice the interactions are screened at 
large distances by remote gates. Instead of modifying the 
Coulomb potential to account for the gate, it will be more 
convenient to simply subtract from Eq. HA1|I a divergent 
constant U{b). Then the series converges, and in the im- 
portant region < 36 the potential can be presented in 
analytic form as 



Uix) = -[F{x/b)-Fil)l 



F{z) 



2i} 



2V' 



3 -f z 



(A2) 
(A3) 



where i}}{z) is the digamma function. 

Evaluation of the energy level splitting AE for a par- 
ticle of mass m in a double- well potential U{x) is a well- 
known problem of quantum mechanics, and the result^ 
is given by 



AE 



TlLJ 



en 



exp 



y/2m[U{x) - huj/2]dx 



(A4) 

Here uj = ^JU"{b)/m is the frequency of small oscilla- 
tions near the minima x = ±b of the potential U{x), and 
X = ±a are the classical turning points at energy Thjj/2, 
i.e., a = b — y^h/muj. 

To evaluate AE with the correct prefactor, one has 
to carefully account for the small energy huj/2 in the 
exponential. The resulting level splitting can be written 
as 

^Vhuj^mb^ e«e-'^°, (A5) 



AE- 



where 



y^2mU {x) dx, 




1 - z 



dz. 



(A6) 



(A7) 
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An alternative solution^ of the problem using the instan- 
ton technique leads to a result that can be also presented 
in the form (|X5|l - (|X7jl . 

In order to apply this result to the evaluation of the 
exchange constant J, one has to keep in mind that x is 
the relative position of two neighboring electrons, x — 
xo — xi, and replace the mass in Eqs. (|A5p . (|A6|) with 
the reduced mass m/2. One then finds Sq — r\\/hl as 
with the numerical coefficient 



y/F{z) - F{l)dz w 2.817. 



(AS) 



Substitution of this result into Eq. l|A5p gives the leading 
exponential behavior of Eq. (|ll|l . 

The numerical parameter ^ defined by Eq. I)A7|I de- 
pends only on the shape of the barrier separating the 
two minima of potential U{x). For the potential ljA2|l . 
we find 




F"{1) 



2[F(z)-F(l)] 1- 



dz « -0.423. 



(A9) 

Substituting this result in Eq. ljA5|l we find the exchange 
constant 



J = 



exp 




(AlO) 



with K « 2.203. Expressing the prefactor in terms of the 
Fermi energy, we obtain Eq. (|12|l . 

It is worth mentioning that because of the singularity 
of the potential U{x) at a; = 0, the validity of the WKB 
approximation used in the derivation^ of formula ljA4|l is 
limited to |x| ^ qb- Moreover, since the potential (|A2p 
is not integrable up to the singularity, it represents an 
impenetrable barrieriS Thus the true value of the level 
splitting for potential ljA2|l is AE = 0. On the other 
hand, the electrons in a quantum wire are not strictly 
one-dimensional due to the finite width w of the wire. 
As a result the singularity of the Coulomb interaction 
potential is cut off at short distances cc ~ w. In GaAs 



devices 



> 



aB, which justifies the WKB approxima- 



tion. In carbon nanotubes it may possible to achieve the 
regime w ^ a^; a more sophisticated approach to the 
calculation of the exchange constant J is required in this 



case 



59 



which shifts the charge field 4>p{x) 4>p{x) + "^^(O- 
As a result the boundary condition (|24|l is replaced with 
(j)p{Q,t) = 0, but the Hamiltonian Q acquires a time- 
dependent perturbation 



V 



ihU^dtU 



q{t) / np(x)dx. 



(B2) 



The perturbation l)B2|l leads to excitation of plasmons 
and to dissipation of energy into the Luttinger liquid. To 
find the energy W dissipated in unit time, it is convenient 
to diagonalize Hp by introducing the plasmon destruction 
operators 



/ e{kx)smkx { -W — ( 
J yTr\j Kp 



bk = I ^ikx) sinkx \ — \l 4'p{x) + i^j -^Ilp{x) j dx, 

(B3) 

where 9{y) is the unit step function. Note that in or- 
der to satisfy the boundary condition 0p(O,i) = the 
wavefunctions of the plasmons were chosen in the form 
ipk(x) = •\/2/7r 0{kx) sin kx] positive and negative k cor- 
respond to excitations to the right and left of the bound- 
ary X = 0, respectively. 

Upon the transformation to the new variables (|B3|I . 
the two terms in the Hamiltonian take the form 



hujkblbkdk, 



V 



ihin COS Lot 



Z^p 



bk 



■dk. 



(B4) 

The perturbation V leads to both emission and absorp- 
tion of plasmons with energy Tiuj. The total energy dis- 
sipated in unit time can be evaluated using the Fermi 
golden rule as 



W 



2lT 



h \2e^2Kp) Tiuo 



[{1 + fk)nLo~ Jkfioj]. (B5) 



Regardless of the values of the plasmon occupation num- 
bers fk, expression HB5|) reduces to W — \l'^Rp with the 
resistance 



APPENDIX C: DISSIPATION OF ENERGY BY A 
SCATTERER IN A FERMI GAS 



APPENDIX B: RESISTANCE OF A QUANTUM 
WIRE 

Let us derive the resistance ()28|l of an infinite quan- 
tum wire in the regime of applied current. The wire is 
modeled by the Hamiltonian Q with the time-dependent 
boundary condition (|24|l . It is convenient to apply to the 
Hamiltonian a unitary transformation 

U = exp ( -i^!^ ( np{x)dx] , (Bl) 



Let us consider the dissipation of energy in a non- 
interacting Fermi gas in the presence of a moving scat- 
terer. We assume that the single-particle Hamiltonian 
has the general form 

H{y,p,t) = Ho{y + qosinujt^p). (CI) 

Here the hermitian operator Hgly^p) is independent of 
the coordinate y in the regions corresponding to the 
leads, y ^ ±oo. The Hamiltonian H35|) obviously sat- 
isfies these conditions for smoothly varying J[y] after the 
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discrete site number I is replaced by a continuous coor- 
dinate y. The y-dependent central part of the Hamilto- 
nian Ho{y,p) can be viewed as a scatterer with energy- 
dependent reflection coefficient Ti{e). Condition IjCip im- 
plies that the position of the scatterer oscillates with am- 
plitude go- 

In the limit of small go one can expand Eq. (jCl|) and 
present the Hamiltonian as 



Here we have approximated the energies near the Fermi 
level as = hvF{\k\ — kp), accounted for the double 
degeneracy of the energy levels Cfc , and introduced 



The matrix element [9y-ffo]fc.fc' is defined as 



(C7) 



H{y,p,t) ^ Ho{y,p) + qo sin ujtdyHo{y,p). (C2) 

The time-dependent perturbation leads to absorption 
and emission of energy quanta huj by the fermions. The 
rates of these processes may be found using the Fermi 
golden rule, and one obtains the increase W of the en- 
ergy of the system in unit time in the form 



W 



27r, 



-huj 



dkdk' 



kk' 



h 

x[/(efc)-/(e'fc)]<5(efc-efc'+M- (C3) 



Here k labels the eigenstates of Hamiltonian Hq with 
energies e^. The occupation numbers of these states are 
given by the Fermi function f{ek)- The eigenfunctions 
have scattering wave asymptotics 



My) = ^ X 

for positive k and 
1 



e'^'^y + TkC-^'^y , at y 



tkc'^y, 



at y 



iky 



tke'^' 



at y 

rke~''^y, at y 



-oo, 

-CO. 



— oo, 

-l-CXD. 



(C5) 



for negative k. Here and tk are the reflection and 
transmission amplitudes; the reflection coefficient is de- 
fined as TZ{ek) = |r-fcp. 

In the limit of low frequency w — > expression (|C3|) 
can be further simplified, 



W 



dek[-nek)][C+(ek) + C-{^k)]. (C6) 



[dyHo]kk'^ J dyrk'iy)[dyHo]My)- (C8) 

Integrating by parts and taking advantage of the fact 
that ipk is an eigenfunction of Hq , we find 



[dyHo]kk' = (cfe -Cfe') J dy^l,{y)dy'ipk{y)- (C9) 

To evaluate C±(£fc) we need to find the divergent at k' — > 
±fc part of the integral in Eq. (|C9|I . Since the divergences 
originate at ?/ ^ ±oo, one can use the scattering wave 
asymptotics (|C4|I and HC5|I in Eq. (jC9|l . This results in 



C+{ek) 



1 



{hvFkFf[n{ek)f 



(CIO) 



C-(efc) = ^{nvFkFf'R{ek)[l~'R{ek)]. (Cll) 



(C4) Substituting these results into Eq. ljC6|l . we find 



W=-{u:qofkl / de[-}'{e)]n{e). 



(C12) 



To apply this result to the evaluation of the spin contri- 
bution to the resistance within the XY model approxi- 
mation, one should substitute ujqo — lo/e and kF = 7r/2. 
Then Eq. takes the form 



W 



R^^ = ^,l de[~ne)\Jl{e). (C13) 



The result for the resistance coincides with Eq. I|4Q|I for 
the appropriate reflection coefficient TZ{e) = 6{\e\ — J). 
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The numerical coefficients in the Hamiltonian @ 
are different from those obtained by conventional 
bosonization— of the Heisenberg model. The discrep- 
ancy is due to the fact that the standard procedure uses 
different bosonic fields (j) = (j)a/V2 and 9 = ^29^. 
One can easily check that the derivation of Eq. Q is valid 
not only for a Wigner crystal, but for any translation- 
invariant model of interacting electrons. 
The asymptotics at j: — > ±cx3 in Eq. 12611 are chosen to 
correspond to outgoing plasmon wave. 
Experiments with cleaved edge overgrowth wiresp show 
quantization below 2e^ /h, but this value is still indepen- 
dent of the electron density. The reduction of the quantized 
value was attributed to the imperfect coupling of the wire 
to the leadsi^ 

It is worth mentioning that the violation of the spin-charge 
separation in Eq. 13011 is due to the inhomogeneity of the 
system. In the homogeneous case J = const, and the spin- 
charge separation is restored. 

The standard solution— applies the WKB approximation 
to the ground states of the oscillatory potentials near the 
minima of U{x). As a result it underestimates the level 
splitting 1A41 by a factor-^ of ^yr/e ^ 1.075. 



